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Abstract 

GRB afterglow is reasonably described by synchrotron emission 
from relativistic blast waves at cosmological distances. We perform 
detailed calculations taking into account the effect of synchrotron self 
absorption. We consider emission from the whole region behind the 
shock front, and use the Blandford McKee self similar solution to 
describe the fluid behind the shock. We calculate the spectra and the 
observed image of a GRB afterglow near the self absorption frequency 
Ua and derive an accurate expression for I'a- We show that the image 
is rather homogeneous for < i^a, as opposed to the bright ring at the 
outer edge and dim center, which appear at higher frequencies. We 
compare the spectra we obtain to radio observations of GRB970508. 
We combine the calculations of the spectra near the self absorption 
frequency with other parts of the spectra and obtain revised estimates 
for the physical parameters of the burst: K 52 = 0.53, Ce = 0.57, es = 
0.0082, ni = 5.3. These estimates are different by up to two orders of 
magnitude than the estimates based on an approximate spectrum. 


1 Introduction 

The detection of delayed x-ray, optical and radio emission following a GRB, 
known as GRB afterglow, is described reasonably well by emission from 
a spherical relativistic shell, decelerating upon collision with an ambient 
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medium (Waxman 1997a, Wijers, Rees & Meszaros 1997, Katz & Piran 1997, 
Sari, Piran & Narayan 1998). A relativistic blast wave expands through the 
ambient medium, continuously heating up fresh matter as it passes through 
the shock. In these models, the GRB afterglow is the result of synchrotron 
emission of the relativistic electrons of the heated matter. 

Several recent works considered emission from various regions on or be¬ 
hind the shock front (Waxman 1997c, Sari 1998, Panaitescu & Meszaros 
1998, Granot, Piran & Sari 1998 (GPS hereafter), Gruzinov & Waxman 
1998). These authors considered the spectra near the peak frequency Vm and 
found that exact calculations of the spectrum could differ by up to one order 
of magnitude from simpler estimates. 

We consider emission from an adiabatic highly relativistic blast wave ex¬ 
panding into a cold and uniform medium. We consider the effect of the whole 
volume behind the shock front, the importance of which was stressed in GPS. 
The hydrodynamics is described by the Blandford McKee (1976 denoted BM 
hereafter) self similar solution. We consider synchrotron emission and we ig¬ 
nore Gompton scattering and electron cooling. Similar to GPS, we consider 
several models for the evolution of the magnetic held. 

Synchrotron self absorption becomes signihcant below a critical frequency 
Va called the self absorption frequency. We assume that I'a ^ ^m, where Um 
is the peak frequency, which is reasonable for the hrst few months after the 
burst. While the spectrum near was quite extensively studied, so far 
only order of magnitude estimates of the spectrum near the self absorption 
frequency were done. 

For a system in which z/^ the spectrum, F^, is proportional to z/^ for 

^ z/a, rather than the standard because almost all the low frequency 
radiation is emitted by electrons with a typical synchrotron frequency much 
higher than z/^ (Katz 1994). The spectra far above z/q is proportional to 
z/^/3_ j ]2 this paper we explore the spectra near z/^, and hnd how the two 
asymptotic forms join together. An analysis of the spectra over a wider range 
of freqnencies was made by Sari, Piran & Narayan (1998), and a detailed 
analysis of the spectra near the peak frequency, taking a full account of the 
BM solution, was made in GPS and Gruzinov & Waxman (1998). 

The physical model is described in §2. In §3 we describe the compn- 
tational formalism. The spectra for several magnetic held models and the 
observed image of a GRB afterglow at various frequencies are presented in 
§4. In §5 we compare the calculated spectra to radio observations of the 
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afterglow of GRB970508. When we use the modihed calculations of the self 
absorption spectra and of the spectra around the peak (GPS) and the cooling 
frequency (Sari, Piran & Narayan 1998) we hud new estimates of the param¬ 
eters of GRB970508. These estimates are different by more than an order 
of magnitude than estimates based on a simpler broken power law spectra 
(Wijers & Galama 1998). 

2 The Physical Model 

The underlying model assumes an ultra-relativistic spherical blast wave ex¬ 
panding into a cold and uniform medium. The blast wave constantly heats 
fresh matter, and the observed afterglow is the result of synchrotron emission 
of the relativistic electrons of the heated matter. We consider an adiabatic 
evolution, where the fluid behind the shock is described by the BM self 
similar solution. It has been numerically verihed that for an adiabatic evolu¬ 
tion with general initial conditions the solution approaches the BM solution 
(Kobayashi, Piran & Sari 1998). The BM solution is expressed in terms of 
the similarity variable y, which is dehned by: 

X = 1 + 167^ , (1) 

where R is the radius of the shock front, r is the distance of a point from the 
center of the burst and 7 / is the Lorentz factor of the matter just behind the 
shock. The BM solution is given by: 

n' = 47/noX"®/^ , 7 = , e' = , (2) 

where n' and e' are the number density and the energy density in the local 
frame, respectively, 7 is the Lorentz factor of the bulk motion of the matter 
behind the shock, rup is the mass of a proton and = rii x lcm“^ is the 
proper number density of the unshocked ambient medium. 

We consider three alternative models for the magnetic held: R, and 
Brad- R satishes e'^, = e^e' (i.e. equipartition) everywhere. For the two 
other magnetic held models we assume they acquire e'^, = e^e' on the shock 
front, and from then on, evolve according to the “frozen held” approximation. 
This implies that the magnetic helds equal where 5 = 7/18 when the 
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magnetic field is in the radial direction {Brad) and 6 = —7/36 when the 
magnetic held is perpendicular to the radial direction {B±) (see GPS). 

We assume that the energy of the electrons is everywhere a constant 
fraction of the internal energy: e'^i = e^e', and that the shock produces a 
power law electron distribution: A^(7e) = for 7e > 7mm- We use 

p = 2.5 wherever a dehnite numerical value of p is needed. The constants 
K and 7 mm in the electron distribution can be calculated from the number 
density and energy density: 

7mm = (- — -r^ ’ K = {p- l)n'7^“/ . (3) 

\p — 1J ri'nieC^ 


We consider frequencies which are much lower than the typical syn¬ 
chrotron frequency u i^syn, where the spectral power of an electron emitting 
synchrotron radiation can be approximated by: 


q^B'sin a / u' 

r(l/3) mec2 \y{yj 



3'jlqeB' sin a 
inrUeC 


(4) 


(Rybicki & Lightman 1979) where P is the gamma function, R' is the mag¬ 
netic held, me and qe are the mass and the electric charge of the electron, 
respectively, and a is the angle between the directions of the electron’s ve¬ 
locity and the magnetic held, in the local frame. 

Ri and 7 ^ are the location of the shock and its Lorentz factor on this 
line of sight and Ri/ji is an estimate of the shell thickness in the local frame. 
The optical depth along the LOS can be approximated by where the 

absorption coefficient a/, is taken at ?/ = y = 1. The “back of the envelope” 
estimate for the self absorption frequency uo is therefore the frequency which 
satishes: al,^Ri/'yi = 1 . uq is given by: 


z/o = 4.24x 10^(1+ ;2)-^ 


3p + 2; 


{p - 1)«/^ 

(p-2) 


1 l/5^1/5 3/5 jt 
-^52 ''-1 ) 


(5) 


where ^ is the cosmological red shift and E = E ^2 x lO^^ergs is the total 
energy of the shell. We also dehne a “standard” flux density Eq, which is an 

^for the energy of the electrons to remain finite we must have p > 2. 
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approximate expression for the flux density at the self absorption frequency. 
The peak flux and peak frequency can be approximated by: 


F 

^ m 


{l + z)niRf 

3 4 WsynJ 


(1/ = X = 1) 


— ^syniSlmini V X 1 ) ) 


( 6 ) 

(see Fq and ut in GPS) where = (4/3)(j'rC7ge^ is the total emitted power 
of an extreme relativistic electron (Rybicki & Lightman 1979), ax is the 
Thomsom cross section and = dx 2 s x 10^®cm is the luminosity distance 0. 
Now we can dehne: Fq = FraiiZa! 


Fq — 1.31(1 + z) 


p + 2' 

3p + 2^ 


1/5 


(p - '-f'X-t ,-i 


(p-2) 


2/5 9/10 7/10^1/2 . 

-^52 J-davs^'^y 


^52 


■ days^ 


(7) 


where T^ays is the observed time in days. 

We express the observed frequency v in units oi vq-. v = 0z/o, thus intro¬ 
ducing the dimensionless variable 0 which we use to express our results. 


3 The Formalism 

We consider a system that is moving relativistically while emitting radiation, 
and obtain a formula for the flux density measured by a distant observer, 
allowing for self absorption. We denote quantities measured in the local 
rest frame of the matter with a prime, while quantities without a prime are 
measured in the observer frame. 

The specihc intensity P (energy per unit time per unit area per unit 
frequency per unit solid angle) satishes the radiative transfer equation: 

^ = F - ajy , (8) 

ds 

where jy and are the emission coefficient and the absorption coefficient, 
respectively, and s is the distance along the beam. When self absorption 
becomes important the contribution to P depends on the optical depth along 
the path within the emitting system and P should be integrated separately 
for every trajectory. 

^Here and latter we use a general cosmological model. For a given model we should 
substitute the appropriate expression for dill + z, Hq). 
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Figure 1: The egg-Shaped curve is the boundary of the region from which 
photons reach a distant observer simultaneously. The observer is located far 
to the right and the symmetry axis is the LOS from the source of the burst to 
the observer. The bold line represents the trajectory of a photon that reaches 
the observer. For a distant observer these trajectories are almost parallel to 
the LOS, and are therefore characterized by their distance R± from the LOS. 
Ji, = 0 at point A, and reaches its hnal value at point B. A photon emitted 
at a point P can be absorbed or cause stimulated emission at any point along 
the trajectory, until it passes the shock front at point B. 


For the BM solution, the radius of the shock front and the Lorentz factor 
of the matter just behind the shock, on the LOS for a given observed time 
T, are given by: 

7; = , Ri = 5.53 x . (9) 

For a given observed time T, we use the coordinates y and x, instead of the 
polar coordinates r and 6*, where y is defined hy y = R/Ri. R is the radius of 
the shock front at the time t at which a photon must be emitted in order to 
reach the observer at the observed time T. A photon emitted at a location 
(r, 6) at a time t in the observer frame, reaches the observer at an observed 
time T given by: 

T = t-—, (10) 

c 

where fi = cos6* (for more details on the coordinates, see GPS). 

Since the observer is distant, the trajectory of a photon that reaches the 
observer is almost parallel to the LOS in the observer frame (see Figure |^). 
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We can therefore parameterize the various trajectories by their distance from 
the LOS. The distance R± of a point from the LOS is given by: 


Rx. = rsme = Riy^l- - XV^ ■ (H) 

The maximal value of R± is obtained for y = 1, ?/ = 5“^/^ and is given by: 

RL,ma. = 3.91 X , (12) 

We express R± in units of R_L^rnax, introducing the dimensionless variable: 


X = 


R± 


55/8 


R±,r. 


y - xy^ 


(13) 


In order to solve equation |^, explicit expressions for the absorption and the 
emission coefficients are needed. We consider an isotropic electron velocity 
distribution, and since P^, ^ oc sin^^^ a we use the averaged value (sin^^^ a) = 
v^r(l/3)/5r(5/6). The total power per unit volume per unit frequency in 
the local frame is given by: 


P' 

W' 





64^13/632/3 (p _ x)5/3 ^V3^^^4/3^8/3^/1/3 

5 r( 5 / 6 ) ( 3 p — l)(p — 2)'^R 2.^29/18-6 


(14) 

Assuming the emission and absorption are isotropic in the local rest frame 
of the matter, j^, = and the absorption coefficient is given by: 



(p + 2) 
STTrUeP^ 



dlePy>,e{le 


N{le) 

7e 


(15) 


8 ^ 7 / 632/3 + 2)(p - 1)8/3 

5r(5/6) (3p- 2) (P- 2)573,5/3^4/3^^,5/3^13/9-5 ’ 

(Rybicki & Lightman 1979). Keeping in mind that j ^/and a^v are Lorentz 
invariant and v' = 1 ^ 7(1 — we can write the radiative transfer equation 
(equation H) explicitly. We hrst write an expression for the optical depth to 
the observer, which will later help us gain some intuition for the results. The 
optical depth to the observer is given by: 



2 rVmaxix) y5/3 

4RR ly ^^yio/9-'5(a;,i/)(l + 7y(a;,i/)i/4)2/3 ’ 
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where xi.^) v) is obtained from equation]^, and ymax{x) is obtained by solving 
Xix,y) = 1 for y. 

We define a “typical” specific intensity Iq by: Iq = S^{y = x = ^), where 
the sonrce fnnction Si, is dehned as: Si, = iyjay. Expressing ly in nnits of 
Jq: ly = -fi/do, we write eqnation in terms of the dimensionless variables Iy\ 

diy / Syx^^^ 

dy _ Txy'^y/^ V(i ~ 



where we took ds = Ridy. This eqnation can be solved nnmerically for a 
given magnetic held model and a given valne of 0, for different valnes of x, 
thns obtaining Iy{x). The observed hnx density is given by: 


Fy{T) ^ 


(1 + ^) 
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I dS±Iy = 27 r(l + z) _ (13) 


As can be seen from equation [T^, the surface brightness is proportional to 
ly-. dFy/dS± oc ly. Thls means that we also obtain the observed image of a 
GRB afterglow, by calculating the observed hnx density in this way. 


4 The Spectra and Observed Image 


Solving equations and 
netic held models {B, Brad 


we calculate the spectra for the three mag- 
and B±). These spectra are shown in Figure 


y. We dehne the self absorption frequency as the frequency at which the 
asymptotic high and low frequency power laws meet. For the equipartition 
magnetic held model B we obtain: 


Ua = 0.247Z/0 = 1.05 X 10^(1+ z) 


-1 


P 


+ 2y'qp-i)s/' 


Z7'^/3„3/5tt 
U/co 77.1 Hz 


3p + 2 


p — 2 


-B 


52 


n 


(19) 

while for the i?_L model it is lower by 6% and for the Brad model it is higher 
by 14%. For n -C the system is optically thick , and therefore the ra¬ 
diation reaching an observer is essentially emitted near the edge facing the 
observer. It rehects the electron distribution there, and it is independent of 
the magnetic held model. For u > Ua the system is optically thin and the 
hnx density is diherent for the various magnetic held models. 













Figure 2: The spectra for different magnetic field models. The frequency Ua 
and the flux density are defined at the point where the extrapolations 

of the power laws meet, as is illustrated for the equipartition magnetic field 
model B. Ua is constant in time, oc and both don’t change 

significantly between the different magnetic field models. For 0 1 (z/ 

Va) the system is optically thick, and the flux density reflects the electron 
distribution (or the Lorentz boosted “effective temperature” of the electrons) 
and is independent of the magnetic field model. 
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Figure 3: The optical depth r to the observer, divided by the maximal 
optical depth (r = Tmax = 1-08 x is black), r^ax is obtained at 

X = Rx!RL,max = 0.93. The contour lines are equally spaced with a 5% 
interval between following contour lines. 


The ratio TjTmax (where r is the optical depth and Tmax is its maximal 
value) for the equipartition B model is shown in Figure |^. Since oc 
everywhere, rjTmax is frequency independent. The contour lines of Ty are 
dense where the absorption coefficient is large. Tmax is obtained at x = 
0.93, i.e. quite close to the edge of the image, since then the whole trajectory 
to the observer is relatively close to the shock front, implying large values of 
ay and a large contribution to Ty. Tmax is a good indicator for the opacity 
of the system, and for the equipartition B model it is given by: Tmax = 
1.08 X {vjvaY^R. 

We dehne Fy^,ext as the extrapolated flux density at z/q (see Figure |^). 
For the equipartition B model we obtain: 


Fu.,ext = O.IO 8 F 0 = 142(1+^) 


p 


1/5 


3p + 2 


( p - l )“'%-2 .-1 

(p-2) 


. 2/5 


^9/10 7/10.^1/2 . 

-^52 ^1 ^daysP'^y i 


( 20 ) 


While for the B^ model it is lower by 12% and for the B^ad model it is 
higher by 23%. The actual flux density at Va is around 35% lower than 
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Fy^,ext- The values of F^^^ext and Ua are useful, since for 

u ua magnitude smaller than Va (assuming z/^ > I'a) the flux density 

is given by: and for z/^ -C z/ -C z/^ it is given by: 

F. = F,^ ,ext{^/■ Both approximations are already good within a few 
percent for frequencies a factor ~ 3 below or above z/q, respectively. 

It is interesting to compare the spectra obtained for the BM solution 
to that obtained for a simplistic model of a static homogeneous disk. This 
comparison should help us learn whether the fact that the spectra is rounded 
up near z/^ should be attributed mainly to the specific hydrodynamics used, 
or whether it is a more general feature of a calculation accounting for self 
absorption. For a static disk we obtain the well known result: 

F. = ( 21 ) 

where the constants F*^ and z/* are determined by the radius and width 
of the disk and by the hydrodynamic parameters of the emitting matter 
within the disk. If these parameters are set so that the two stared quantities 
equal those in equations |2^ and |^, respectively, the resulting flux density 
is very similar to that obtained for the BM solution. Substituting equations 
|2y and [T^ into equation ^ can thus serve as a good approximation (better 
than 3%) for the observed flux density at z/ -C Vm- This similarity implies 
that the shape of the spectrum near Va is rather independent of the specific 
hydrodynamic solution considered. On the other hand, the exact values of 
Va and F^^^ext depend signihcantly on the hydrodynamics, and can not be 
determined without the detailed calculation. A simplihed calculation for u* 
could for example yield uq (equation |^) instead of Ua that is given in equation 

m 

The surface brightness as a function oi x = Rx/Rx,max for the equipar- 
tition magnetic held model B is shown in Figure for a few representative 
values of 0 = zz/z/q. An illustration of the observed image of a GRB after¬ 
glow, which is implied from the surface brightness, is shown in Figure]^. As a 
reference, the image for v S> zz^, which was derived in GPS, is also presented. 

For u Ua the system is optically thin and the result coincides with that 
obtained in GPS for zz -c zz^. There is a bright ring near the outer edge of 
the image, and the surface brightness at the center of the image is 58% of its 
average value. As zz decreases, the contrast between the edge and the center 
of the image decrees and the system becomes increasingly optically thick. 
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Figure 4: The surface brightness divided by the average surface brightness, 
as a function of x, for Lo5fio(0) = —1.5, —1 — 0.75, —0.5, 0 (corresponding to 
v/ua = 0.13,0.40,0.72,0.28,4.05). At high frequencies the contrast between 
the center and the edge of the image is larger than at low frequencies. At low 
frequencies the system is optically thick, and the surface brightness reflects 
the “effective temperature” of the electrons at the edge of “egg” depicted in 
Figure |^, on the side facing the observer. 
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Figure 5: The observed image of a GRB afterglow, at several frequencies. 
At high frequencies there is a bright ring near the outer edge of the image, 
and the contrast between the center and the edge of the image is larger than 
at low frequencies. At low frequencies the surface brightness increases as 
one moves from the center towards the edge, until it drops very sharply, due 
to the fact that the system becomes optically thin near the edge. The last 
image, for v 3> is taken from GPS and it is brought as a reference, to 
illustrate the change in the relative surface brightness along the image over 
a large range of frequencies. At z/ 3> there is a thin bright ring at the 
outer edge of the image, and the surface brightness at the center is only a 
few percent of its maximal value. 
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As a result, the surface brightness reflects the Lorentz boosted “effective 
temperature” of the electron distribution, at the outer edge of the “egg” 
depicted in Figures and |^, on the side facing the observer. Larger values 
of X correspond to smaller shock radii R, and since the electrons posses a 
larger Lorentz factor (i.e. a larger “effective temperature”) at smaller radii, 
the surface brightness increases with x. This is true as long as long as the 
optical depth is still large. Since the length of the trajectory within the 
system approaches zero as x —> 1, for every given frequency u the system 
becomes optically thin for values of x sufficiently close to 1. For frequencies 
smaller than by more than one or two orders of magnitude, this occurs at 
1 — a; -C 1 so that the drop in the surface brightness near a: = 1 is extremely 
sharp. For z/ -C the surface brightness at the center of the image is 77% 
of its average value, resulting in an almost uniform disk, rather than a ring 
which is obtained for u > z/^. The uniformity of the image for u Ua stands 
out even more when compared to z/ S> z/^, where there is a thin bright ring 
on the outer edge of the image and the surface brightness at the center is 
only a few percent of its maximal value. 

5 Comparison to Observations 

We now fit the theoretical spectra calculated, to the radio frequencies obser¬ 
vations of the afterglow of GRB970508 (Fig. 4 of Shepherd 1998). Perhaps a 
future near enough GRB will result in a sufficient resolution, that will enable 
us compare the predicted images of a GRB afterglow to observations. 

Since both z/q and Fq depend on the parameters of the model (see equa¬ 
tions ^ and 0), we have two degrees of freedom in trying to fit the calculated 
spectra to the observed data. We fit the equipartition magnetic field model 
B to the data in Shepherd et al (1998), and obtain: 

z/, = 3.1±0.4x lO^Hz , F,^,ext = 450 ± 37/iJy , (22) 

with y^/dof = 0.48. The fit is presented in Figure ^ 

The errors quoted are statistical errors. The actual errors are probably 
larger due to uncertainty in the radio flux due to scintillation, and due to 
the fact that not all the observations are exactly simultaneous (see Shepherd 
1998). 
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Figure 6: A fit of the calculated spectra, to radio observations of the afterglow 
of GRB970508. We obtained: Va = 3.1 ± 0.4 x lO^Hz and = 450 ± 

37/iJy. 


Substituting these results in equations and ^ respectively, we obtain 
two constraints on the physical parameters of the burst. The luminosity 
distance depends on the cosmological model. For hi = 1, A = 0 (which 
implies: di = 2(1 + z — \/l + z)c/Hq), Hq = 65Kms“^Mpc“^ and z = 0.835 
we find: 


(p - ifF / p + 2 
(p-2) \3p + 2 

(p - ifF / p + 2 \ 
(p-2) [ 3 P + 2 ) 



1 l/ 5 pl /5 3/5 
-^52 ^1 


12/5 p9/10 
-^52 


n 


7/10 

1 


5.43 , 

1.22 . 


(23) 

(24) 


Substitution of the value for Va from the data (equation |^) into equation 
|T^, yields an equation different by a factor of ~ 2 from equation 22 of Wijers 
& Galama (1998). The difference is mainly due to a different theoretical 
value they used for They also used a slightly different value of Va as 
corresponding to the same observational data. This factor of ~ 2 implies 
signihcant corrections to the values of the physical parameters of the burst 
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TABLE 1 


Estimates of the Physical Parameters of GRB970508. 


model 

broken 
power law 

modihed 

Va 

modihed z/q, 

I Fum 1 

E52 

3.7 

2 

0.53 

ee 

0.13 

0.24 

0.57 

es 

0.068 

0.011 

0.0082 

ni 

0.035 

0.70 

5.3 


Table 1: The first row depicts the values of the physical parameters of 
GRG970508 as calculated by Wijers & Galama (1998). The other rows show 
how these values change when some of the equations they used in the cal¬ 
culation are corrected. The hrst line lists the measurable quantities, whose 
equations were altered. 

(see the second row of Table 1). For example, ni (or n in their notation) 
becomes a factor of ~ 20 larger: ni = 0.7 instead of ni = 0.035. 

We modify the estimates further using more detailed calculations of the 
spectrum near the peak flux (GPS), and use a different estimate for the 
cooling frequency Vc (Sari, Piran & Narayan 1998). We discover that the 
values obtained for the physical parameters of a burst are very sensitive to 
the theoretical model of the spectrum. In order to illustrate this we show 
in Table 1 the physical parameters of GRB970508 for different estimates of 
the spectrum. The hrst column uses a broken power law spectrum (Wijers 
& Galama 1998); the second column uses a corrected theoretical value for Va 
from equation and a corrected observational value for from equation ^ 
the third column adds a modihed theoretical values to Vm and (which are 
taken from GPS - denoted there as Ppeak and Fy^max) and a different estimate 
for the cooling frequency (from Sari, Piran & Narayan 1998), keeping the 
observational values from Wijers & Galama (1998). 

Our best values (with all modihcations added) are E 52 = 0.53, Cg = 0.57, 
es = 0.0082 and rii = 5.3. These values differ by one to two orders of 
magnitude than the values obtained by Wijers & Galama (1998) using a 
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broken power law spectra. One must keep in mind that it is difficult to 
obtain an accurate estimate of the various observables from the observed 
data. It is especially difficult to determine and Vc- Therefore, more 

than obtaining a better estimate of the physical parameters of the burst, 
these calculations show the sensitivity of this method and the need for more 
accurate data. 

Sufficient observational data has been gathered on the radio afterglow 
of GRB980329, to enable a £t similar to the one we made for GRB970508. 
Such a fit was carried out by Taylor et al (1998). The theoretical formula for 
the flux density that was used for the £t is identical to equation 1^. Since 


this is a good approximation for the shape of the spectra near z/q, we can 
use the values extracted from the data to obtain constraints on the physical 
parameters of GRB980329. 

The values used by Taylor et al for the flux density are mean values over 
the first month. Since Fy^ oc this corresponds to the flux density at 

~ 15 days. The values extracted from the data are: 


Va = 1-3 X 10 


10 


Hz 


F,.„ext = 600/xJy . 


(25) 


Substituting these results in equations and respectively, we obtain 
two constraints on the parameters of GRB980329. For G = 1, A = 0 and 
Hq = 65Kms“^Mpc“^ we hnd: 


2 (p - ifF / p + 2 ■ 

(1 + ;^) (p-2) 1^3p + 2 


3/5 


1 l/5pl/5 3/5 
G -^52 ^1 


i" ( P±li 


e-B E. 


^ + z — 1J (p — 2) y 3p + 2 

where the red shift 2 : of this burst is not yet known. 


52 


26 , (26) 
1.5 , (27) 


6 Discussion 

We have considered synchrotron emission from a system moving relativisti- 
cally, taking into account the effect of synchrotron self absorption. We have 
assumed an adiabatic evolution and used the self similar solution of Bland- 
ford Sz McKee (1976) to describe the matter behind the shock. This solution 
describes an extreme relativistic spherical blast wave expanding into a cold 
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uniform medium. Our calculations accounted for the emission from the whole 
region behind the shock front. 

We have assumed a power law distribution of electrons with an isotropic 
velocity distribution. Three alternative models have been considered for the 
evolution of the magnetic held, including an equipartition model. We have 
calculated the hux density at frequencies near the self absorption frequency 
Vai under the assumption that i^a ^ ^m- 

We have obtained an expression for the self absorption frequency Ua, 
which we dehned as the frequency at the point where the extrapolations 
of the asymptotic power laws above and below the Ua meet. The value we 
obtained for Ua is close to the value obtained by Wijers & Galama (1998), 
and a factor of ~ 4 for p = 2.5 (or a factor of ~ 7 for p = 2.2) larger than 
the value obtained by Waxman (1997b). 

We have calculated the observed spectra near for three different mag¬ 
netic held models (Figure The spectra dihers more than a few percent 
from the asymptotic power laws only for frequencies less than half an order 
of magnitude above or below z/^. This applies to all the magnetic held models 
we have considered. 

Together with the afterglow images obtained in GPS we now have a com¬ 
plete set of the observed image over a wide range of frequencies (see Figures 
^ and 1^). For u ^ Um there is a thin bright ring on the outer edge of the 
image, and the surface brightness at the center of the image is only a few 
percent of its maximal value. For Ua ^ ^ there is a wider ring and the 
contrast between the center and the edge of the image is smaller, with 58% 
of the average surface brightness at the center. For z/ -C z/^ the image is even 
more homogeneous, with 77% of the average surface brightness at the center. 

The observed image in the radio frequencies is of special importance, since 
the best resolution is obtained in radio, with VLBI. A sufficient resolution 
could be reached with a nearby GRB {z ~ 0.2) to resolve the inner structure 
of the image in radio frequencies. 

We have htted the theoretical spectra to observational data of the af¬ 
terglow of GRB970508 (Shepherd 1998) and extracted the values of z/^ and 
Fua,ext from the data (equation |2^). Substituting these values into the the¬ 
oretical expressions we obtained two equations for the physical parameters. 
The behavior of the spectra near z/^ can supply two other equations, namely 
the equations for the peak frequency and for the peak flux (adding only one 
independent equation to the two equations for z/q and Fy^^ext)- An additional 


18 


independent equation is obtained for The power law p of the electron 

distribution can be determined from the high energy slope {v S> Vm)- It 
is therefore possible with sufficient observational data to determine all the 
physical parameters of the afterglow: p, eg, e^, E ^2 and rii. A similar cal¬ 
culation (using an equation for the cooling frequency Pc instead of 
was made by Wijers & Galama (1998) for GRB970508. In fact combining 
this equation we could even have an over-constrained system and check for 
consistency of the solution. The equation we obtained from differs by a 
factor of ~ 2 from the corresponding equation given in Wijers & Galama. 
This has a significant effect (up to a factor of ~ 20) on the values of the phys¬ 
ical parameters of the burst which they have given. We have shown (Table 
1) that if other equations are corrected as well, the values of the physical 
parameters vary even more (up to two orders of magnitude). 
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